> 

o 



'Td 



X 



M- furcations in coupled maps 

Sang-Yoon Kim* 

Department of Physics 

Kangwon National University 

Chunchon, Kangwon-Do 200-701, Korea 

Abstract 

We study the scaling behavior of M-furcation (M = 2,3,4,...) sequences of 



^ . M"-period (n = 1,2,...) orbits in two coupled one-dimensional (ID) maps. 

<^ . 

^^ . Using a renormalization method, how the scaling behavior depends on M 

is particularly investigated in the zero-coupling case in which the two ID 
0> . maps become uncoupled. The zero-coupling fixed map of the M-furcation 

^.^1 renormalization transformation is found to have three relevant eigenvalues 5, 



a, and M (6 and a are the parameter and orbital scaling factors of ID maps. 



o 

,^ . respectively). Here the second and third ones, a and M, called the "coupling 

O 

eigenvalues", govern the scaling behavior associated with coupling, while the 

first one 6 governs the scaling behavior of the nonlinearity parameter like the 

case of ID maps. The renormalization results are also confirmed by a direct 

numerical method. 
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I. INTRODUCTION 

Universal scaling behaviors of M-furcation (M = 2,3,4,...) sequences of M"-cycles 
(n = 1,2, . . .) (i.e., M"-period orbits) have been found in a one-parameter family of one- 
dimensional (ID) unimodal maps with a quadratic maximum. As an example, consider the 
logistic map 

xt+i = f{xt) = 1- Ax'^t, (1) 

where t denotes the discrete time. As the nonlinearity parameter A is increased from 0, a 
stable fixed point undergoes the cascade of period-doubling bifurcations accumulating at a 
finite parameter value Aoo{= 1.401 155 . . .). The bifurcation sequence corresponding to the 
MSS (Metropolis, Stein, and Stein [Q) sequence i?*" (for details of the MSS sequences and 
the (*)-composition rule, see Refs. [0,|1) exhibits an asymptotic scaling behavior 0. 

What happens beyond the bifurcation accumulation point A^q is interesting from the 
viewpoint of chaos. The parameter interval between A^o and the final boundary-crisis point 
Ac{= 2) beyond which no periodic or chaotic attractors can be found within the unimodality 
interval is called the "chaotic" regime. Within this region, besides the bifurcation sequence, 
there are many other sequences of periodic orbits exhibiting their own scaling behaviors. In 
particular, every primary pattern P (that cannot be decomposed using the (*)-operation) 
leads to an MSS sequence P*"". For example, P = RL leads to a trifurcation sequence 
of 3"-cycles, P = RL^ to a tetrafurcation sequence of 4"-cycles, and the three different 
P = RLR^, RL^R, and RL^ to three different period-5" sequences. Thus there exist 
infinitely many higher M-furcation (M = 3,4,...) sequences inside the chaotic regime. 
Unlike the bifurcation sequence, stability regions of periodic orbits in the higher M-furcation 
sequences are not adjacent on the parameter axis, because they are born by their own tangent 
bifurcations. The asymptotic scaling behaviors of these (disconnected) higher M-furcation 
sequences characterized by the parameter and orbital scaling factors, 6 and a, vary depending 
on the primary pattern P MR 



In this paper we consider two symmetrically coupled ID maps. This coupled map may 
help us to understand how coupled oscillators, such as Josephson-junction arrays or chem- 
ically reacting cells, exhibit various dynamical behaviors |T2|-p^. We are interested in the 
scaling behaviors of M-furcations (M = 2, 3, . . .) in the two coupled ID maps. The bifurca- 
tion case (M = 2) was previously studied in Refs. |IB|-pU[. Here we extend the results for 



the bifurcation case to all the other higher mult if ur cat ion cases with M = 3,4,... in the 
zero-coupling case where the two ID maps become uncoupled. In Sec. Jl|we investigate the 
dependence of the scaling behavior on M using a renormalization method. It is found that 
the zero-coupling fixed point of the M-furcation renormalization transformation has three 
relevant eigenvalues S, a, and M. The scaling behavior associated with coupling is governed 
by two coupling eigenvalues (CE's) a and M, while the scaling behavior of the nonlinearity 
parameter is also governed by the eigenvalue 6 like the case of ID maps. As an example, we 
numerically study the scaling behavior associated with coupling in the trifurcation sequence 
in Sec. |T^ and confirm the renormalization results. Finally, a summary is given in Sec. ^. 

II. RENORMALIZATION ANALYSIS 

In this section we first introduce two coupled ID maps and discuss stability of orbits, 
and then study the scaling behavior of M-furcations (M = 2, 3, . . .) in the zero-coupling case 
using the renormalization method developed in Refs. |[T5| , p^ . It is found that there exist 



three relevant eigenvalues 6, a, and M. As in the case of ID maps, the scaling behavior 
of the nonlinearity parameter is governed by the eigenvalue 6, irrepectively of coupling. 
However, the scaling behavior associated with coupling depends on the nature of coupling. 
In a linear-coupling case, in which the coupling function has a leading linear term, it is 
governed by two CE's a and M, whereas it is governed by only one CE M in the other cases 
of nonlinear coupling with leading nonlinear terms. 

Consider a map T consisting of two symmetrically coupled ID maps. 



yt+1 = F{yt,xt) = fiyt) +g{yt,xt), 

where f{x) is a ID uniinodal map with a quadratic maximum at x = 0, and g{x,y) is a 
couphng function. The uncoupled ID map / satisfies a normahzation condition 

/(O) = 1, (3) 

and the couphng function g obeys a condition 

g{x, x) = for any x. (4) 

The two-coupled map (^ is invariant under the exchange of coordinates such that x ^^ y. 
The set of all points which are invariant under the exchange of coordinates forms a symmetry 
line y = X. An orbit is called an "in-phase" orbit if it lies on the symmetry line, i.e., it satisfies 

Xt = yt for all t. (5) 

Otherwise it is called an "out-of-phase" orbit. Here we study only in-phase orbits, which 
can be easily found from the uncoupled ID map, Xt+i = f{xt), because of the condition (H). 
Stability of an in-phase orbit with period p is determined from the Jacobian matrix J of 
T^, which is the p-product of the Jacobian matrix DT of T along the orbit: 



f[DT{xt,Xt) = f[ 
t=i t=i 



^ f\xt) - G{xt) G{xt) ^ 



(6) 



\^ G{xt) f\xt) - G{xt) ) 

where the prime denotes a derivative, and G{x) = dg{x,y)/dy \y=x', hereafter, G{x) will be 
referred to as the "reduced coupling function" of g{x,y). The eigenvalues of J, called the 
stability multipliers of the orbit, are: 

Ai = n/'(^t)' X2 = f[[f{xt)-2Gixt)]. (7) 

t=i t=i 

Note that the first stability multiplier Ai is just that of the uncoupled ID map and the 

coupling affects only the second stability multiplier A2. 



An in-phase orbit is stable only when the moduli of both multipliers are less than or 
equal to unity, i.e., — 1 < Aj < 1 for i = 1,2. A tangent (period-doubling) bifurcation occurs 
when each stability multiplier A, increases (decreases) through 1 (—1). Hence the stable 
region of the in-phase orbit in the parameter plane is bounded by four bifurcation lines 
associated with tangent and period-doubling bifurcations (i.e., those curves determined by 
the equations Aj = ±1 for i = 1,2). 

We now consider the M-furcation renormalization transformation TV, which is composed 
of the M-times iterating (T^*^-*) and rescaling (B) operators: 

Af{T) = BT^^'^B-\ (8) 

Here the rescaling operator B is: 



B 



la ,\ 



V a I 



(9) 



because we consider only in-phase orbits. 

Applying the renormalization operator TV to the coupled map @) n times, we obtain the 
n-times renormalized map T„ of the form, 

. 1 ^*+i = Pn{xt,yt) = fnixt)+gnixt,yt), 
yt+1 = Fn{yt,xt) = fniyt) + gniyt,xt). 

Here /„ and Qn are the uncoupled and coupling parts of the n-times renormalized function 
Fn, respectively. They satisfy the following recurrence equations: 

fn-,i{x)= a/W(-), (11) 

a 

a a a 



gn^i{x,y)=aF^^^\^,l)-ae^\^), (12) 



where /W(x) = Uli''-'Kx)) and F(^\x,y) = F„{F(^-'\x,y), F(''-'\y,x)). 

The recurrence relations ([TT|) and (|12|) define a renormalization operator TZ of transform- 
ing a pair of functions (/, g) : 



(f \ 

Jn+1 
\ 9n+l J 



n 



Jn 

\gn) 



(13) 



The renormalization transformation TZ obviously has a fixed point (/*, g*) with g*{x, y) = 0, 
which satisfies 7l{f*,0) = (/*,0). Here /* is just the ID fixed function satisfying 



nx) = ar(''^{-) 



,x. 



a 



(14) 



where a = l//*(^^ ^''(l); due to the normahzation condition (|]), and it has the form, 



f*{x) = 1 + c^x" + c^x^ + ■ 



(15) 



where c*'s {i = 1,2, . . .) are some constants. The fixed point (/*, 0) governs the critical be- 
havior near the zero-coupling critical point because the coupling fixed function is identically 
zero, i.e., g*{x,y) = 0. Here we restrict our attention to this zero-coupling case. 

Consider an infinitesimal perturbation [h, ip) to the zero-coupling fixed point (/*, 0). We 
then examine the evolution of a pair of functions (/* -|- h, ip) under TZ. Linearizing TZ at 
the zero-coupling fixed point, we obtain a linearized operator C of transforming a pair of 
perturbations {h,(f)): 



h. 



n+l 



\ (.\ ( 



c 



y V^n+l j 



iifi 



\'^n J 
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Cr 
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where 



(16) 



in+l{x) = [Cuhn]{x) 



a5/f^)(^)^a[/W(^)-/:W(^)]Hnear 



a 



a 



a 



a a a 



ipn+i{x,y) = [CcVn]ix,y) 



«^[i^r (-, ^) - e'\-)] - a [Fr)(-, ^) - /rH-)]nnea. 



X,, r„/n,T\,X y 

a a' ' " a" ' " ' a a 



a' 



«r'(/ 



a 



-))'5[i^r 



, /r-^^(-)] 

a a a 



*{M-i)(^\ f*{M-i)(y_ 



:-)j 



a 



-))• 



a 



(17) 
(18) 
(19) 
(20) 
(21) 

(22) 



Here the variations Sfj^^\-) and (5[F^*^^(-, ^) — /^^'*(-)] are introduced as the hnear terms 
(denoted by [/W(f) - f<'^'\^)U.r and [FW(f,f) " /^Hf )]iinear) in h and ^ of the 
deviations of fi^\^) and F^^^)(f , ^) - /i^'^Hf ) from f<^\l) and 0, respectively. A pair 
of perturbations (/i*, </)*) is then called an eigenperturbation with eigenvalue z/, if it satisfies 



(lA 



V 



C 



V^V 



(23) 



V^V 



I.e., 



h*{x) = [Cuh*]{x), 
^*{x,y) = [Cc^*]{x,y). 



(24) 
(25) 



The eigenperturbations of the linear operator C can be divided into two classes. The first 
class of eigenperturbations are of the form [h*, 0) . Here h* (x) is an eigenfunction of the linear 
"uncoupled operator" £„ satisfying Eq. (p^ , which is just the eigenvalue equation in the 
uncoupled ID case. It has been found in Refs. [|,||,|I that there exist a unique eigenfunction 
h*{x) with (noncoordinate change) relevant eigenvalue 6, associated with scaling of the 
nonlinearity parameter. 

The second class of eigenperturbations have the form {0,ip*), where ^p*{x) is an eigen- 
function of the linear "coupling operator" Cc satisfying Eq. (|25|). However, it is not easy 
to directly solve the coupling eigenvalue equation (pSj). We therefore introduce a tractable 
recurrence equation for a "reduced coupling eigenfunction " of ip*{x,y), defined by 

d<^*{x,y) 



<^*(x) = 



dy 



(26) 



y=x 



Differentiating Eq. (^) with respect to y and setting y = x,we obtain an eigenvalue equation 
for a reduced linear coupling operator C^- 



I/$*(X) = [£e$*](x) 



^i^i"^(-) = [i^i"^(-)]nnear 



a 



a 



a a' 



*{M-iy(^^^*fMM-i)'X 



a 



«))• 



(27) 
(28) 

(29) 



Here F{x,y) = f*{x) + (p*{x^y), F2 (x) is a "reduced function" of F^^^^^x^y) defined by 
F2 (x) = dF^^'^\x^y)/dy\y=x., and the variation 5F2 (-) is also introduced as the hnear 
term (denoted by [Fi^'^^(f )]iinear) in $* of the deviation of Ff ^^(f ) from 0. 
In the case M = 2, the variation 5F2 {-) of Eq. ( PH]) becomes 



5Ff (-) = $*(-)/*'(/*(-)) + r'(^)r(r(-)). (30) 

a a a a a 



Substituting 5F2 (-) into Eq. (^), we have 6F2 (-) for M = 3, which consists of three 
terms, 



5Ff (-) = $*(^)/*'(r(^))/*'(r(2)(i)) 

a a a a 

+r'{-W{r{-))r\r^'\-)) 

a a a 



+r(-)r'(r(-))$*(r^'n-))- (31) 

a a a 

Repeating this procedure sucessively, we obtain 5F2 (-) for a general M, composed of M 
terms, 

M-l 



= r(-)r(^^-i)'(r (-)) + ••• 



a a 



where f^^\x) = x. 

Using the fact that /*'(0) = 0, it can be easily shown that when x = 0, the reduced 
coupling eigenvalue equation (^9]) becomes 



M-l 

A$*(0) = [n/*'(/*^*nO))]$*(0). (33) 



Differentiating the ID fixed-point equation (|1^) with respect to x and then letting x ^ 0, 
we also have 

M-i f*'(T^ 

nnr«(0)) = lim^ = «. (34) 



Then Eq. (|33| ) reduces to 

A$*(0) =a$*(0). (35) 

There are two cases. If the couphng eigenfunction ip*{x,y) has a leading hnear term, its 
reduced couphng eigenfunction $*(x) becomes nonzero at x = 0. For this case $*(0) 7^ 0, 
we have the first CE 

ui = a. (36) 

The eigenfunction ^l{x) with CE ui is of the form, 

$*(x) = l + a^x + a^x^H , (37) 

where a*'s (i = 1, 2, . . .) are some constants. For the other case $*(0) = 0, it is found that 
f*'{x) is an eigenfunction for the reduced couphng eigenvalue equation (p9|). Since Eq.(p2D 
for the case $*(x) = f*'{x) becomes 

5Fi^)(-) = MrW(-), (38) 

a a 

the reduced coupling eigenvalue equation ( P^D reduces to 

uf'ix) = Mf'ix). (39) 

We therefore have the second relevant CE 

z/2 = M, (40) 

with reduced coupling eigenfunction $2(^) = f*'{^)- It is also found that there exists an 
infinite number of additional (coordinate change) reduced eigenf unctions f*'{x) [f*\x) — x'] 
with irrelevant CE's a"' (/ = 1,2, . . .), which are associated with coordinate changes. We 
conjecture that together with the two (noncoordinate change) relevant CE's (z/i = a, 1^2 = 
M), they give the whole spectrum of the reduced linear coupling operator Cc of Eq. ( JI7\ ) 
and the spectrum is complete. 
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In order to examine the effect of CE's on the stabihty multiphers of periodic orbits in the 
M-furcation sequences, we consider an infinitesimal couphng perturbation g{x, y) = eip{x, y) 
to a critical map at the zero-coupling critical point, in which case the two-coupled map has 
the form, 

T : < ^°= (41) 

[ yt+i = F{yt,xt) = fA(M){yt) + g{yt,xt), , 

where A^-* denotes the accumulation value of the parameter A for the M-furcation case, 
and e is an infinitesimal coupling parameter. The map T at e = is just the zero-coupling 
critical map consisting of two uncoupled ID critical maps. It is attracted to the zero-coupling 
fixed map consisting of two uncoupled ID fixed maps under iterations of the M-furcation 
renormalization transformation A/" of Eq.(|^). 

The reduced coupling function G{x) of g{x, y) is given by [see Eq. 

dip{x,y) 



Q{x) = e^{x) = e- 



dy 



(42) 



y=x 

The nth image $„ of $ under the reduced linear coupling operator Cc of Eq. (|27|) is of form, 

$„(a;) = [C:^]ix) 

~ aiz/"$^(a;) + a2^'2/*'(a;) for large n, (43) 

because the irrelevant part of $„ becomes negligibly small for large n. Here ai and a2 are 
some constants. 

The stability multipliers Ai^„ and A2,n of the M^-cycle of the map T of Eq. (|4T| ) are the 
same as those of the fixed point of the n-times renormalized map A/'"(T) [0, which are 
given by 

Al,n = fLi^n), A2,n = fLi^n) - 2G„(x„). (44) 

Here /„ is the uncoupled part of the nth image of {fAM),g) under the renormalization 
transformation TZ, G„(x) is the reduced coupling function of the coupling part gn{x,y) of 
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eiz/" + e2Z/^ for large n. (46) 



the nth image, and x„ is just the fixed point of fn{x) [i.e., £„ = /„(£„)] and converges to 
the fixed point x* of the ID fixed map f*{x) as n ^ oo. In the critical case (e = 0), A2,n 
is equal to Ai,„ and they converge to the ID critical stability multiplier A* = f*'{x*). Since 
Gn{x) ~ [£"G'](a;) = £$„(a;) for infinitesimally small e, A2,n has the form 

A2,n ^ Ai,„ - 2£$„ 

~ A* + £ [eiz/" + 622/2] for large n, (45) 

where Ci = —2ai^l{x*) and 62 = — 2a2/*'(a;*)- Hence the slope Sn of A2,n at the zero- 
coupling point {e = 0) is 

e=0 

Here the coefficients Ci and 62 depend on the initial reduced function $(a:), because the 
constants ai and 0:2 are determined only by $(a;). Note that the magnitude of slope Sn 
increases with n, unless both coefficients ei and 62 are zero. 

We choose monomials a;' (/ = 0, 1, 2, . . .) as initial reduced functions $(a;), because any 
smooth function $(a;) can be represented as a linear combination of monomials by a Taylor 
series. Expressing $(a;) = x' as a linear combination of eigenfunctions of Cc, we have 

<l>(x) = x' = ai<l>i(a;) + a2f*\x) 

+j:M*'{x)[r\x)-xl (47) 

1=1 

where cti is nonzero only for / = 0, and hence zero for / > 1, and all jSi^s are irrelevant 

components. Therefore the slope Sn for large n becomes 

I eia" + esM" for / = 0, 
^n ^ (48) 

I e2M" for / > 1. 

There are two kinds of coupling. In the case of a linear coupling, in which the coupling 
function ip{x,y) has a leading linear term, its reduced coupling function $(a;) has a lead- 
ing constant term. However, for any other nonlinear-coupling case, in which the coupling 
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function has a leading nonlinear term, its reduced coupling function contains no constant 
term. Hence it follows from Eq. (^) that the growth of Sn for large n is governed by the 
two relevant CEs z^i = a and z/2 = M for the linear-coupling case (1 = 0), but by only the 
second relevant CE 1/2 = M for the other nonlinear-coupling cases (/ > 1). 

III. NUMERICAL ANALYSIS 

Taking the trifurcation case with M = 3 as an example, we numerically study the scaling 
behavior associated with coupling in the two coupled ID maps ( Pl) with f{x) = 1 — Ax^ 
and ^{x, y) = — (y™ — x™) {m = 1, 2, . . .), and confirm the renormalization results. In this 
trifurcation case, we follow the S'^-cycles up to level n = 9 and obtain the slopes of Eq. ( ^61) at 
the zero-coupling critical point (A(^\ 0) (A^^^ = 1.786 440 255 563 639 354534447...) when 
the reduced coupling function $(x) is a monomial x' (/ = 0, 1, . . .). 

The renormalization result implies that the growth of the slopes Sn is governed by one 
CE z/2 = 3 for the nonlinear-coupling cases with / > 1, i.e., the sequence of Sn obeys a 
one-term scaling law asymptotically: 

Sn = dirl (49) 

where di is some constant. We therefore define the growth rate of the slopes as follows: 

ri,n = -^- (50) 

Then it will converge to a constant ri asn -^ 00. As an example consider the case ^(x) = x. 
Figure |l| shows three plots of A2,n(^^'',^) versus e for n = 5,6, and 7. For e = 0, A2,n 
converges to the ID critical stability multiplier A* (= 1.872 705 929 . . .) as n ^ 00. However, 
when e is nonzero it diverges as ra — * 00, i.e., its slope Sn at the zero-coupling critical point 
diverges as n ^ 00. The sequence {ri „} of the growth rate of Sn is shown in the second 
column of Table |. Note that it converges fast tori = 1/2 = 3. We have also studied two other 
nonlinear-coupling cases with / = 2, 3 and found that the sequences of ri „ also converge fast 
to ri = 1^2 = 3. 
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However, in a linear-coupling case with / = 0, two relevant CE's ui = a {= —9.277341 . . .) 
and z/2 = 3 govern the growth of the slopes Sn- We therefore extend the simple one-term 
scaling law (HOf) to a two-term scaling law: 



Sn = c?ir" -|- d2r2 for large n, (51) 

where di and d2 are some constants, and |ri| > |r2|. This is a kind of multiple scaling law 



21|,|22|. The equation (plD gives 

Sn+2 = QlSn+l — ?2>S'n, (52) 

where qi = ri + r2 and g2 = ''"i'"2- Then ri and r2 are solutions of the following quadratic 
equation, 

r^ - qir + q2 = 0. (53) 

To evaluate ri and r2, we first obtain gi and g2 from S^s using Eq. (^2[): 

Q Q C C C2 C Q 

_ *J7i+l*Jn '~>n+2i~>n-l _ '~>n+l '-'n'Jn+2 (KA\ 

yi Q2 Q C ' ^^ "c2 C C ■ \ ) 

Note that Eqs. (pT|)-(|5^) are valid for large n. In fact, the values of g^'s and r^'s (i = 1,2) 
depend on the level n. Thus we denote the values of g^'s in Eq. (U) explicitly by gi,„_i's, 
and the values of r^'s obtained from Eq. (|53D are also denoted by r^ „_i's. Then each of them 
converges to a constant as n ^ oo: 

lim qi^n = qi, hm Vi^n = n, i = 1, 2. (55) 

When $(x) = 1, plots of A2,n(^^\^) versus e for n = 2,3, and 4 are shown in Fig. |^. 
The slopes S^ at (y4|^\ 0) obeys well the two-term scaling law (|5T|). Sequences {ri^„} and 
{^2,n} are shown in the third and fourth columns of Table |. Note that they converge fast 
to ri = ui = a and r2 = 1^2 = 3, respectively. 
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IV. SUMMARY 

The scaling behavior of M-furcations is studied in two symmetrically coupled ID maps. 
Using a renormalization method, the dependence of the scaling behavior on M is particularly 
investigated in the zero-coupling case. It is found that the zero-coupling fixed map of the 
M-furcation renormalization operator has three relevant eigenvalues 5, a, and M . As in the 
case of ID maps, the eigenvalue 5 governs the scaling behavior of the nonlinearity parameter, 
irrespectively of coupling. However, the scaling behavior associated with coupling depends 
on the nature of coupling. In a linear-coupling case, it is governed by two CE's a and M, 
whereas it is governed by only one CE M in the case of a nonlinear-coupling case. Taking 
the trifurcation case as an example, we also study the coupling effect on the second stability 
multipliers of 3"-cycles by a direct numerical method and confirm the renormalization results. 
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TABLES 
TABLE L In a nonlinear-coupling case $(x) = x, a sequence {ri^„} for a one-term scaling law 

is shown in the second column, and in the linear coupling case ^(x) = 1, two sequences {ri^„} and 

{'^2,n} for a two-term scaling law are shown in the third and fourth columns, respectively. 

<l>(x) = X <l>(x) = 1 

n n^n n^n r2,n 

1 2.997 929154 -9.276 54316 2.005 8 

2 3.000141141 -9.277 415 78 2.692 5 

3 2.999 990 417 -9.277 335 54 2.927 8 

4 3.000 000 651 -9.277 34150 2.984 5 

5 2.999 999 956 -9.277 34109 2.996 7 

6 3.000 000 003 -9.277 34112 2.999 3 

7 3.000 000 000 

8 3.000 000 000 
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FIGURES 

(3) 

FIG. 1. Plots of the second stability multipliers A2,n(^oo , e) versus e for n = 5,6,7 in a non- 
linear-coupling case <I>(a:) = x. 

FIG. 2. Plots of the second stability multipliers A2,n(^oo , e) versus e for n = 2,3,4 in the 
linear-coupling case <I>(x) = 1. 



